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Abstract. 

Conformally flat spacetimes with an elastic stress energy tensor given by a diag- 
onal trace- free anisotropic pressure tensor are investigated using 1+3 formalism. 
We show how the null tetrad Ricci components are related to the pressure com- 
ponents and energy density. The 1+3 Bianchi and Jacobi identities and Einstein 
field equations are written for this particular case. In general the commutators 
must be considered since they supply potentially new information on higher or- 
der derivatives of the 1+3 quantities. We solve the system for the non rotating 
case which consist of ODEs of a spatial coordinate. 

1 Introduction 

The theory of elasticity in the context of general relativity was developed in the 
mid twentieth century. The need for such a theory came in the late 1950s with 
Webers bar antenna for gravitational waves [T], in order to explain how these 
waves interact with elastic solids. Actually, for this phenomena the weak-field 
approximation was sufficient in the treatment of the problem given by Weber. 
Only in 1973, in a paper by Carter and Quintana [2], did a fully developed non- 
linear theory of elasticity adapted to general relativity appear, remaining to this 
day a standard reference in the field, although the basic theoretical framework of 
this theory had already been given by Souriau in [3J. Also before the article by 
Carter and Quintana, work by Maugin made considerable contributions to the 
field @]-[5]. Lately, the theory of elasticity in general relativity was reconsidered 
by Magli and Kijowski [6]- [7] and Chritodoulou jS], in this work they explore 
the gauge character of relasticity. Authors such as Beig and Schmidt, have 
proven several existence and uniqueness theorems [9] . More recently, Karlovini 
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and Samuelsson have given a self contained formulation of general relativistic 
elasticity in [10] .They applied the theory of elasticity to spherically symmetric 
space-times and studied radial and axial perturbations [H]-[I3]- Park [T3] estab- 
lished existence theorems for spherically symmetric static solutions for elastic 
bodies and Brito, Varot and Vaz [T5] have obtained static shear-free and non- 
static shear-free solutions for spherically symmetric elastic spacetimes. Calogero 
and Heinzle |16j studied the dynamics of Bianchi type I elastic spacetimes. 

In this work we study a very simple elastic model with a diagonal trace-free 
anisotropy pressure tensor using the 1+3 extended frame approach and follow- 
ing the notation given in Uggla [T7]. We note that work on the covariant 1+3 
splitting of fluid spacetime geometries was first initiated by Eisenhart and Synge 
and continued by Godel, Raychaudhuri and other authors such as Schiicking, 
Ehlers, Sachs and Triimper (related references [18], [19], [20], [21]). In the paper 
by Uggla the basic dynamical equations of the extended 1+3 orthonormal frame 
approach are explicitly given in terms of variables that are naturally adapted 
to the 1+3 structure, and they include the Bianchi and Jacobi identities, the 
Einstein field equations and the commutators. This formulation is analogous 
to the Newmann Penrose approach [22] in the sense that a null congruence is 
replaced by a timelike congruence. The general properties of the 1+3 orthonor- 
mal frame can be found in books such as Wald [53] and Felice and Clarke [2"4] 
and in Edgar [55] . Several applications have been discussed in Pirani [55] , Ellis 
[57] and MacCallum [55]. A more complete list of references can be found in 

The organization of this paper is as follows. In section 2 we outline the theory 
of the 1+3 formalism and present the 1+3 split of the commutators, curvature 
variables and their field equations, namely Bianchi and Jacobi identities and 
Einstein filed equations. We will follow the same notation convention for tensor 
indices used in [17j . Spacetime coordinate tensor indices will be denoted by 
letters from the second half of the Greek alphabet (/it, v, p, ... = — 3) while 
spatial coordinate indices are represented by letters from the second half of 
the Latin alphabet (i,j,k,... — 1,...,3). Orthonormal frame indices will be 
denoted by letters from the first half of the Latin alphabet (a, b, c, ... = 0, 3) 
while spatial frame indices are chosen from the first half of the Greek alphabet 
(a,/3,7,... = 1 — 3). In section 3 we relate the Ricci spinor components to 
the 1+3 quantities. In section 4 we study conformally flat spacetimes with 
an elastic source given by a diagonal trace-free anistropy pressure tensor using 
the 1+3 formalism. We determine the 1+3 equations for a very simple case of 
non rotating solutions. In particular we study the non rotating case where the 
system obtained is an ODE system of a spatial coordinate. 

2 1+3 Formalism 

When studding a dynamical model in general relativity possessing an energy- 
momentum-stress tensor with a timelike eigendirection, for example perfect flu- 
ids, the associated timelike vector field u on the spacetime (M,g) determines 
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the projection tensors U and h, which project parallel and orthogonal to u in 
the tangent space at each point p e (M,g), respectively, u is chosen to be a 
unit timelike vector: 

«/ = -!. (!) 
and the projection tensors U and h are defined by 

U» v = -u»u v (2) 



fc" v = 5"„ + (3) 

Due to the existence of this singled out timelike direction u a covariant 1+3 
tensor decomposition of all geometrical objects of physical value can be made 
using U and h. 

We will denote covariant derivative by V and the totally antisymmetric 
permutation tensor by tj^'" 7 . 

The well known kincmatical fields associated with the timelike congruence 
u are defined by 

ii^ = h^ v u p V p u v (4) 



e = h^v^u" (5) 



[^(^%) - 3 V^KVv) (6) 



u lu , = -h p ^h <r v] (V p u a ), (7) 

where w/ 1 is the acceleration vector, the rate of expansion scalar, cr^is the 
rate of shear tensor and uj^ u is the vorticity tensor. The magnitude of the rate 
of shear a, the vorticity vector w M and the magnitude of the vorticity oj are 
defined as 

* 2 - \^^ v (8) 

^ = ^rr^uvpu* (9) 
1 



u? = -lo^lo^ = W/i oA (10) 



The vector field u is hypersurface forming if u> = 0. 

In the orthonormal frame approach one chooses at each point of the space- 
time manifold (M, g) a set of four linearly independent 1-forms e a such that the 
line element is given by 
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ds 2 = rj ab e a e b , (11) 

where "q a b — diag[— 1, 1, 1, 1] is the constant Minkowskian frame metric. The 
vectors e a represent the dual basis. 

In the 1+3 orthonormal frame formalism one aligns the timelike direction 
of the orthonormal frame with the tangent of the preferred timelike congruence 
(e = u). 

2.1 Commutators 

The commutator functions are defined by 

[e a ,e b ] = l c ab e c , (12) 

where the frame vectors e a act as differential operators on geometrical objects. 
The Ricci rotation coefficients T a b c can be written in terms of the commutation 
functions as 

7 a 6c = r a ch - r\ c . (13) 

The commutator functions with one or two indices equal to zero can be ex- 
pressed in terms of the frame components of the kinematic quantities associated 
with the timelike congruence. The angular velocity Q a is defined by 

n a = lv abcd eb-Vu(e c )u d . (14) 

The purely spatial components of the commutation functions 7 c * ( g 7 can be 
decomposed |29j a a and n a p by means of 

1 a pi = 2a [0 8 a l] +e MS n Sa , (15) 

where n a p is symmetric and e^ 7 is the totally antisymmetric three-dimensional 
permutation tensor. 

The commutators are given by expressions for 7 a 6c and their 1+3 decompo- 
sition results in the following equations 

[e 0) ej = u a e -[^®8f i a +cj f) a + e f) ai {u"t -W)]efi (16) 
[e a ,ep] = -2e a ^ 7 w 7 e + pa^fJ 7 ^ + e a/3 5n 57 ]e 7 . (17) 

2.2 Curvature variables and their field equations 

The relation between the Riemann curvature tensor and the Ricci rotation co- 
efficients is given by 

R a bcd = e c (r a i,d) — ed(r a bc) + r a ec r e fcrf — r a e dT e b C — Y a i >e v e C d 1 (18) 
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and the 16 Jacobi identities are then given by the following equations 

R a [bcd] = 0; Rabcd = Rcdab- (19) 

The Einstein field equations with a non zero cosmological constant can be 
written as 

Rab = T ab - ^Tg ab + Ag ab , (20) 

and the 1+3 decomposition of the energy- momentum-stress tensor T ab with 
respect to the timelike vector field u in the fluid description of phenomena of a 
matter source is given by 

Tab = fJ,u a u b + 2q {a u b) + ph ab + n ab , (21) 

where \i denotes the total energy density scalar, p denotes the isotropic pres- 
sure scalar, q a denotes the energy current density vector, and ir ab denotes the 
anisotropic pressure tensor. Note that 

q a u a = 0; TT ab u b = Q; ir a a = 0; n ab = ir ba . (22) 

The matter fields need to satisfy an appropriate thermodynamical equation 
of state in order to describe the physics of the underlying fluid spacetime geom- 
etry under consideration. 

The completely tracefree Weyl conformal curvature tensor is defined by 

C ab cd = R ab cd - 2S^ [c R b ^ d] + \m a [c 5 b d] . (23) 

When there is a singled out timelike direction u, it is convenient to decom- 
pose the Weyl tensor into its electric part given by 

E ab - C cedf h c a u e h d b u f , (24) 

and magnetic part 

H ab - *C cedf h c a u e h d b u f , (25) 
where *C abcd is the dual of C abcd defined by 

1 , 

*C abcd = -rj ab e: 'C e f cd . (26) 

The electric and magnetic tensors are symmetric and tracefree and satisfy 
E ab u b = H ab u b = and lead to the expression 

C ab cd = [^ a e 5 b ^S^ [c S h d] - V ab efV o h cd }E e g ^u h 

- 2[ v ab ef 5\ c 5 h d] +8^ e 5 b ^ f ^ h cd ]H%ufu h . (27) 
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Considering equation ([23]) and using ([20)1 , (|2"T1) and ([27)1 we obtain the fol- 
lowing expression for the Riemann curvature tensor 



R ab cd = [46KS b] f 5\ c 5 h q-if b ef if h cd \E e g uJu h 
- 2[r, ab ef S^ c S h d] +5^ e 6 b ^rf h cd ]H e g ufu h 
+ 25^ [C [^+P)u b] u d] + g% + u% + A,] 

+ ^ + A)S a [c S b d] . (28) 

Inserting expression ([TBI for the Riemann tensor into the left hand side of 
l]28p one obtains the 10 Einstein field equations, the 16 Jacobi identities and 
expressions for E a b and H a b in terms of the basis variables 

{e, a a , n a p, fi, p, q a , TT a p}, 

and their eo and e a frame derivatives 

Einstein field equations 

e (e) = ~^6 2 + (e« +u a - 2a a )(u a ) - 2a 2 + 2uj 2 

- I( M + 3p) + A (29) 
e (a afi ) = -Q<7 a ' 3 + {8~< {a e~ l +u {a + a {a ){u^) + 2^ a tt^ + tt q/3 

- *S a0 - ^ Q/3 [(e 7 + u 7 + a 7 )(?i 7 ) + 2w 7 fF] 

+ e^ a [20 7 <7« < 5-n^ (5 u 7 ] (30) 
= /i- ie 2 + o- 2 - w 2 - 2uj a n a - l*R + A (31) 

= (e^-3a^(a^)-^^e^(e) + n Q ^ + g Q 



e a ^[(ep + 2up - a/OK) + 7^ 5 a 5 7 ], (32) 



where 



*S' Q/ 3 = e( a (a /S )) + 6 Q/ 3 - -S a p[e 1 (a 7 ) + 6 7 7 ] 

- £ 7<5 (a(e| 7 | - 2a hl )(n p)s ) (33) 

*ii = 2(2e a -3a a )(a a )~^b a a (34) 

6 Q(3 = 2n ai n 1 p - n? \n a p. (35) 

Jacobi identities 
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e (a a ) = -l(5 afi ep + u a +a a )(Q) + ^(ep + up~2ap){a afi ) 

- ^e a ^(ep+u f3 -2a fj )(^-n j ) (36) 

- 2n fj \(uj s - n s )} (37) 
eoK) - ~Gu a + a a pu^ + \n a pv! 5 - e a ^[^(ep - a fj )(u 7 ) 

+ upSly] (38) 

= (ep - 2a p )(n a P) - ^Qu a - 2a a + e a ^[ep(a 7 ) 

+ 2ujpVL 1 ] (39) 

= (e a - u a - 2a a )(uj a ). (40) 

The Bianchi identities are given by 



V [a i? bc]de = -V [a [R c b] - ^6 c b] R] = -V [a [T c b] - ±6 c b] T\, (41) 

and can be expressed in terms of E ab and H ab and the 1+3 basis variables 
associated with T ab 

Bianchi identities 



e (M) = -(Li+p)6-(e a + 2u a -2a a )(q a )-a a pTr a P (42) 

e (q a ) = -^eq a -S^e (p)-^+p)u a 

- (e p +up- 3a )(n^) - o a pq? + e a ^[(cop 

] (43) 



6 

+ I„7 7 ff^ + ij^[i( e7 + 2ti 7 + a 7 )(gT) 
- 3cr 7(5 (E lS -\^ 5 ) + Zn. (S H~< s ] + e^ a 

[(e 7 + 2?i 7 - a 7 )(H^ s ) - K - 2fi 7 )(£7« 4 
+ + - 3n^H^ (44) 
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- ^ af> ) + 3n (a 7 (£^ - \^ Y <) - 6 a P[a l5 H^ s 
\^<P + n, s (E^ s - - e^«[(e 7 - a 7 ) 

+ {u 7 -2Q 7 )H^ 5 } (45) 
= {ep - 3ap)(E a ^ + - l -5^e^) + ±Gq a 

- \° a m P + 3u H a P - e a ^[a f3S H\ + |^g 7 

+ n ps {E s 1+ l -^ 1 )\ (46) 
= (ep-3a fj )(H a P)-(i2 + p)Lo a -3cup(E a P 

+ apsiE^ + ^^-npsHW. (47) 

3 Relation between the covariant 1+3 and null 
tetrad decomposition of the curvature tensor 

3.1 Relasticity 

The formulation of the theory of relativistic elasticity is based on a configuration 
mapping tp : M — > X from the spacetime M to the three-dimensional material 
space X, which represents the collection of particles of the material and is 
equipped with a material Riemannian metric 7. The material coordinates will 
be denoted by y A , A = 1, 2, 3, and let represent the coordinates in M. This 
mapping describes the configuration of the material and gives rise to a rank three 
matrix (jj A ^) , p € M, y A ^ — which is called the relativistic deformation 
gradient. The velocity field of matter u M a future oriented, timelike unit vector 
field , which spans the one-dimensional Kernel of the relativistic deformation 
gradient is defined by the conditions y^u^ = 0, u^u^ = — 1. The pulled-back 
material metric 

V = V 1AB = V*V*1/AB, (48) 

can be used to construct the relativistic strain tensor 

where h^ v = + u^u Vl which measures the state of strain of the material. If 
s^v = 0, the material is said to be in an unstrained state. 
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k^ v has three positive eigenvalues n\, ri|, n§. The particle density of the ma- 
terial is defined by n = \fkF\, — n\n 2 n^ and Ui, 71.2, 71.3 are called linear particle 
densities. The energy density can be defined by p = ne, where e represents 
the energy per particle. A constitutive equation for the material is given by 
specifying p as function of three scalar invariants of k^ v , p = p(Ii, h, h), if 
we assume that the internal energy of an elastic deformation accumulated in 
an infinitesimal portion of the material is invariant with respect to the space- 
time orientation of the material. For the scalar invariants one can choose 
h = h = h = Kk%k\. Since n 2 = \ (/j 5 - 3hh + 2/ 3 ) , 

one can use the particle density n as one of the scalar invariants. 

The stress-energy tensor for elastic matter is given by 

T M „ = pu^u v + p^, (50) 
where the pressure tensor p^ v is defined by 

p^=php V +n l w, (51) 

p being the isotropic pressure and Tr a b the anisotropic pressure tensor. 

An equation of state that is compatible with the stress energy tensor (150|) 
can be written as 

p = p + pa 2 , (52) 

where p is the unsheared energy density p is the modulus of rigidity, a 2 the 
shearing scalar and 

p = p+(Q-l)a, P = n 2 ^, n = lf, (53) 

an p an 

being e the energy per particle. 

3.2 Determining the NP curvature components in 1+3 

If we represent the orthonormal tetrad as {e a M }, the relation between the NP 
null tetrad vectors fc M , ra M , fh 11 } and the 1+3 vectors are given by 



J" = ^(eo^-ea") 

= -^(^W) 

= l = (e^+ie 2 ' 1 ), (54) 

and its dual 
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1 

1 

-y=(e ft — ie fj) 

-j=(e\ + ie 2 p). (55) 



fc P = ^(e^ + e 3 ^) 



The 4- velocity is equal to eo M = u M = ^(fc^ + an d the metric will then be 
given in terms of this basis by 

= 2m (Al m M) - 2l^k^. (56) 



Contracting the Einsteins field equations (|20l) in coordinate components with 
the appropriate null vectors and considering the momentum tensor (|21|1 in coor- 
dinate components we obtain the relation between the Ricci spinor components 
and the 1+3 quantities 



$00' = $22' = |(M+P + ^33) (57) 
$11' = ^(M+P + ^n +^22 -^33) (58) 
$01' =$10' = |(-7ri3+wr 23 ) (59) 

$12' =$21' = |(7Tl3-Mr 2 3) ( 60 ) 
$02' = $20' = ^(TTll - 7T22 ~ «2tT 12 ) (61) 

A = ^(M-3p). (62) 

On the lefthand side we have the Ricci spinor components and on the righthand 
side we have the 1+3 components of the symmetric anisotropic pressure tensor 
7r a f,, the energy density scalar /i and isotropic pressure scalar p. If we consider 
7r a b to have zero components then we have a perfect fluid source and the relation 
between the Ricci spinor components and the 1+3 quantities are those obtained 
in [3D]- 

Defining the tensor 

Qfj,u — C 'pjy paU h ' v? — E^p + iHpp, (63) 
the dual part of the Weyl tensor can be expressed in terms of Qp V as 
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*C pvpa u v u a — 8u[ fl Q ll ][ (T u p ]+2g fi y ) Q a j v — 2g u y,Q^ fM 

+ 2ie f _ l „~ fS u 1 U[ p Q cr ] S + 2ie pai su 1 u^Q u ^ . (64) 

Considering that the spinor equivalent of *C pvpa is tA'B'^CD'^ abcd and 
contracting (f64|) with the null tetrad vectors one can express the Weyl spinor 
components in terms of the 1+3 components of Q pv as 



*0 = 


5«u- 


Q 2 2™ 




*1 = 




-Qu) 




*2 = 








*3 = 




-Q13) 




*4 = 




Q22 + 





4 Conformally flat metrics with elastic stress 
energy tensor 

It is well known that the Weyl tensor C pvpG vanishes iff the spacetime is con- 
formally flat. By definition, the metric of a conformally flat spacetime can be 
written as: 

ds 2 = F 2 {t,x,y, z)(-dt 2 + dx 2 + dy 2 + dz 2 ). (66) 

All conformally fiat solutions with a perfect fluid, an electromagnetic field, 
or a pure radiation field are known. 

Given that the Ricci spinor components transform under null rotations in 
the following way 



$00' 


-> $00' 








$01' 


-> $01' - 


+ a$oo' 






$02' 


-> $02' 


■f 2a$oi' 


+ a 2 $oo' 




$11' 


— S- $n' 


■f a$io' -+ 


- a$ i' 4 


aa$oo' 


$12' 


-> $12' 


■f 2a$ii' 


+ a 2 $i ' 


+ a$02' + 2aa$oi' + a 2 a$oo- 


$22' 


-> $22' 


■f 2a$2i' 


+ 2a$i2< 


+ 4aa$ii' + 2aa 2 $oi' 




+ 2a 2 a 


$10' + a 2 


$20' + oe 


2 $ 02 ' +a 2 a 2 $ 00 ', 



with a a complex function, if 7Ti3, 7r23, 7Ti2 ^ and 7r 2 3 = ^(n + p + tt33)(itii — 

27T 2 2 - /i-p), vr 2 3 = |(/i+p + 7T33)(-27ril +7T 2 2 -fl-p), 7T 2 2 = i (tTh - 27T 2 2 - 
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u — p)(— 7Tn + 27T22 — M — p)j one can choose a basis (a = — such that in 
that basis 7Ti3, 7r23, 7Ti2 = 0. 

Here we study the special case of conformally flat spacetimes with elastic 
energy-momentum tensor (|50[) . with an anisotropic traceless pressure tensor 
TTjui, that can be put in the form 

7TU = 7T22 = - 2 ^33 : (68) 

with all other components being null and with zero curvature scalar 

A= ^( M -3p)=0, (69) 

and therefore the only non-zero curvature terms are 

$00' = ^O-TTn) = $22' (70) 

$ii' = ^(p + tth). (71) 

If we calculate the Bianchi identities for this special case we obtain the 
following equations for the 1+3 directional derivatives of p and ir as well as 
conditions involving p, 7rn, n a @, a a , U a , uj a , a a/ 3, Q a , 







eo(p) 


4 

= --p6 - 2(Tii7rii 


(72) 






eo(TTll) 


= -4auP + {cru - -6)7rn 


(73) 






ex(p) 


= -(01-^23)^11 


(74) 






ei(7r u ) 


= (ai-n 2 3)7Tii 


(75) 






e 2 (p) 


= -(02 +ni3)7TH 


(76) 






e2(7Tll) 


= (02 + "13)7111 


(77) 






e3(p) 


4 2 

= -g«3P+ gW37Tll 


(78) 






e3(7Tn) 


4 , 2 . 
= 2 U3P + ( 3 ° 3 ~ 3 M 3)ttii 


(79) 


<iuip + (iii - 


- 3ai 4 


- 3n 23 )7Tii 


= 


(80) 


4ii2P + (ii2 - 


- 3a 2 - 


- 3ni 3 )7Tii 


= 


(81) 


2<Ji3P+ ^(3w 2 


-6Q 2 


+ Cl3)7Tll 


= 


(82) 


1 . 

2cr 2 3P - t( 3w 1 


-6fii 


- 023)111 


= 


(83) 


—4u!ip + 


2>i- 


3CT23)7ril 


= 


(84) 
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-4w 2 p + ^(w 2 + 3cri 3 )7Tii = (85) 

^3 = "33 = "12 = 0"12 = (86) 

Oil = 0"22 = _ 2 CT33 

"n = "22- (88) 

We take eri 3 = cr 23 = ^>i = ^2 = ^3 — ^1 = ^2 = ui = ^2 = and 
"23 = Oi; "13 = — a-2 so that (|80"j) - (f8"5)) are automatically satisfied and the 
Bianchi identities reduce to 



eo(p) - -3^e-2crii7rii (89) 

eo(7rii) = -4cri 1 jj+(o-ii--e)7rn (90) 

ei(p) = (91) 

ei(7r n ) = (92) 

e 2 (p) = (93) 

ea(7r u ) = (94) 
4 2 

e 3 (p) = --u 3 p+ -u 3 7Tii (95) 

4 2 

e3(7Tn) = o^3P+ (3a 3 - -u 3 )7Tii. (96) 



As for the Ricci and Einstein field equations they are given by 



e (a 3 ) = -^u 3 @ - ~a 3 Q - u 3 au - a 3 aii (97) 
3 1 3 

ei(ai) + e 2 (a 2 ) + e 3 (a 3 ) = -p - -9 2 + -g\ x + 2a\ + 2aj 

2 6 2 

+ ~a§ (98) 
e 3 (a 3 ) - e (<Tu) - -e 3 (u 3 ) = 8crn - 7Tn + i-£t§ 

+ ^03«s+P- g© a + oii + ol (99) 

e 3 (oii) + ^e 3 (9) = 3a 3 crii (100) 

e 3 (^ 3 ) + e (n n ) = -u 3 fi 3 + 2annn - ^Qn n (101) 

e (e)-e 3 (u 3 ) = -^e 2 - 6o? 1 -3p + u|-2a 3 u 3 (102) 

ei("n) - 2e 3 (a 2 ) = 2a\nu - 2a 3 a 2 (103) 

2e 3 (ai) - e 2 (nn) = 2a 2 rin + 2a 3 ai (104) 



13 



eo(ai) - 7:62(^3) 



-ai8 - oi<tii - a 2 3 
-026 — aaaii + 01^3 



(105) 



eo(a 2 ) + 2 e i(^s) 
ei(u 3 ) = e 2 (u 3 ) 
ei(a 3 ) = e 2 (a 3 ) 

ei(o-n) = e 2 (crn) 
ei (e) = e 2 (6) 



0. 















(106) 

(107) 
(108) 
(109) 
(110) 



In the following section we study the ODE solutions within this particular 
case. 

4.1 Non rotating ODE solutions 

We will start by investigating the case where eo, &\, e 2 acting on all quantities 
are zero in which case the previous equations will reduce to ODEs. A simple 
study of the Bianchi identities shows that eo, e±, e 2 acting on p and ttu are zero 
implies that O = — 3o"n or irn = —Ap. However, if we apply commutator (|16l) 
with a = 3 to all quantities involved p, ttu, an, nu, ai, a 2 , 03, SI3, <d the 
conclusion is that & — 6(Jn otherwise all such quantities are constant also (|17|) 
with a — 1 and /3 = 2 implies that n\\ =0 , (|17[) with a — 1 and /? = 3 implies 
that ai = and (jTTJ) with a = 2 and = 2 implies that a 2 = otherwise all 
quantities are constant. Therefore two cases arise. 

Case AI 7rn + 4p = 9 — 6un = a\ = a 2 = nu = 

The Jacobi identities and the Einstein equations in summary give the fol- 
lowing ODE system 




(111) 



(112) 

(113) 
(114) 
(115) 
(116) 



u 3 = -a-3 



63(0-11) = a 3 CTn 
e 3 (a 3 ) = -9<7? x +2a§ 



If we now choose the coordinates (t, x, y, z) and the 1+3 basis 




(F(z), 0,0,0) 
(0,F(^,0,0) 
(0,0,F(z),0) 
(0,0,0, F(z)), 



(117) 
(118) 
(119) 
(120) 
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its solution for an, a 3 , and F in this basis can be written as 



<Jii = u u (z) (121) 

a 3 = ±{Aa 2 11 +9)i<j 11 (122) 

dz 

fl 3 = Bel^ dz dz, (124) 

where A and B are constants and an is a non constant function of z. If in 
particular the solution is shearless and expansionless, equation (|115[) gives the 
following solution for 03 



fl3 = jrpTTc' (125) 

where C is a constant. 

We note that these equations do not allow a perfect fluid solution since 
7TU = -4p = 0. 

Case A2 an = = «n = 04 = 122 = 

In this case the solution has zero expansion, is shearless and the only extra 
condition given by the Jacobi identities is in relation to the non zero component 
of the angular velocity 



e 3 (ft 3 ) = -^3- (126) 

The Einstein field equations give a relation for the non zero anisotropic 
pressure component, the isotropic pressure scalar, non zero component of accel- 
eration and the non zero component of the geometric object a a 

1 1 2 . . 

Tin = -p- -a 3 + a 3 u 3 , (127) 

along with the non linear ODE system in p, 7Tn, 03, u 3 



1 2 

&zip) = -u 3 p - -u 3 a 3 + -a 3 ul (128) 

63(^3) = 3p~ul + 2a 3 u 3 (129) 

3 3 
— V H — ( 

r 2 



e 3 (a 3 ) = 7;P+n4- (130) 



We note that if 03 = then the equations above give p = tth = and 
therefore this solution reduces to the case of a dust solution and the only ODE 
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one must solve relates to U3. In the case of zero acceleration U3 = then the 
isotropic pressure p is zero and the system reduces to the ODE in a 3 . Also if 
we put TTn in the equations above we obtain p = so therefore this case does 
not allow for a perfect fluid. 

If we choose to solve this system using the basis (|117[) - (|120p the relations 
between the Ricci spinor components and the 1+3 quantities give 7Tn = — 4p 
since $00' = —5^11' so that the equations then give W3 = —03 or 713 = ^03. 
When U3 = — 03 the solution can be obtained by putting a\\ = in equations 
(|llip - (|116[) and is given by (|125[) . The solution when u 3 = ^03 is given by 

" = S-^U + D- (131) 
with D a constant. The solution for ft^ is the same as in (|124j) . 



5 Future Work 

A more complicated problem would be to consider a non-rotating, non accel- 
erated solution with null angular velocity (ui a = u a = Cl a = 0), purely spatial 
components of the commutator functions given by 



n af) = (132) 
a a = (0,0,03), (133) 

and diagonal shear 










(J 11 








-2(711 





The Bianchi identities in this case become 



(134) 



4 

e o(p) = -gP©- 2CTH7TH (135) 

eofau) = -4cthp+ (cr n - ^©)7Tn (136) 

e3(7rn) = 3a 3 7Tii 

ei(p) = e 2 (p) = e 3 (p) = 

ei(7rn) = e 2 (7ru) = 0, (137) 

and the Ricci and Einstein field equations 
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eo(a 3 ) 
e 3 (a 3 ) 
eo(o-n) 



63(0-11) + ge 3 (9) 

e (e) 

ei(a 3 ) = e 2 (a 3 ) 
ei(cru) = e 2 (crii) 
ei(0) = e 2 (e) 

ei(n 3 ) = e 2 (n 3 ) = e 3 (n 3 ) 

One must apply the commutators to the quantities p, 7rn, cth, f2 3 , a 3 , 9. 
In particular the only commutator not identically satisfied when applied to these 
quantities is (|16p with a = 3. This commutator can in fact give rise to further 
information on their second order covariant derivatives. If we write the previous 
equations in a basis such as (|117j) - (|120|) we see that we are dealing with a non 
linear PDE system which of course is very difficult to solve so that we should 
consider generalizing the algorithm described in the ODE case to the study of 
this more general case. 

In the study of non conformally flat sapcetimes with elastic stress-energy 
tensor it would be very useful to have a 1+3 version of the Karlhede classification 
of metrics since there are several techniques for obtaining exact solutions of 
Einsteins equations involving tetrad formulation of these equations in which the 
Karlhede classification is a useful instrument. 
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-ga 3 9 - a 3 an 



— a 3 

2 6 



1 3,1, 
7;P + ttii + 2°"n +2 a3 " 

0CT11 

3a 3 CTn 

_i 2_ 6CTll 2_ 3p 




0. 
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